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1 Introduction 



In the search for a supergravity model dual to a "realistic" strong coupling gauge dynamics (for 
a review see |lj), one should be able to incorporate quarks in the fundamental representation of 
an SU(N C ) gauge theory. Most of the known supergravity backgrounds duals of confining four 
dimensional gauge theories either do not incorporate quarks at all or admit quarks in the adjoint 
rather than the fundamental representation 1 . 

Since the early days of strings it has been understood that fundamental quarks should cor- 
respond to open strings. In the modern era of closed string theory this obviously calls for D 
branes. Certain basic objects of gauge theories like baryons [2J, instantons, monopoles, domain 
walls and others jl] were shown to correspond to wrapped D brane probes. It is thus natural 
to wonder, whether one can consistently add D brane probes to supergravity backgrounds duals 
of confining gauge theories, which will play the role of fundamental quarks. In case that Nf 
the number of D brane probes is much smaller than N c , one can convincingly argue that the 
backreaction of the probe on the bulk geometry is negligible. 

It is well known that open strings between parallel Nf D7 and N c D3 play the role of flavored 
quarks in the SU(N C ) gauge theory on the D3 4d world volume gauge theory. Karch and Katz 
[3] proposed to elevate this brane configuration into a supergravity background by introducing a 
D7 brane probe into the AdS$ x S 5 background. 

This idea was further explored in jB]. In [7] the spectrum of an M = 2 SU(N C ) SYM with 
fundamental hypermultiplet was extracted from the supergravity background of AdS$ x S 5 with 
D7 brane probe . It was found that for massive quarks the spectrum was discrete with a mass 
gap. They also analyzed the semiclassical rotating open strings attached to the D7 brane and 
discussed meson meson interaction. For a related work, see also jH] 

The purpose of this paper is two folded: (i) to find a supergravity dual of a confining gauge 
theory with flavors and (ii) to compute the mass spectra of the pseudo and vector mesons made 
of the flavored quarks from the SUGRA background. The first goal was in fact proposed but 
not explicitly proven in |S] and as for the second goal recall that the computation in [2j involves 
M = 2 SQCD which is not a confining gauge theory. The main idea behind our work is to 
introduce D7 brane probes into the Klebanov Strassler (KS) |S] background (this idea was first 
given in 5J). We will see that this configuration yields an M = 1 SU(N C ) gauge theory with 
massive Nf flavors. The spectrum of the pseudo scalar and vector mesons can be extracted from 

1 Even bifundamental quarks in models of SU(N) x SU(N) turn into an adjoint ( and a singlet) once the 
symmetry group is broken to the diagonal SU (N) group. 
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the spectrum of scalar and vector fluctuations of the Nf D7. It turns out that it is advantageous 
not to use the formulation of the deformed conifold of JU] but rather the coordinates introduced 
in [H]. In the latter picture there is a separation of the coordinates of the three manifold and 
the S 2 . It would be quite interesting to compare our results with a field theory analysis. We will 
leave it as a future problem. 

To understand the construction of the probe D7 branes recall that in the type IIA T dual 
picture the singular locus of the conifold takes the form of a pair of perpendicular NS five branes 
[T2*l I13j. Once the conifold is deformed there will be a "diamond" structure at the intersection 
of the NS5 branes j!4j . The fractional D3 branes of the KS model are mapped by the T duality 
to D4 branes connecting the two NS branes. Now in this type IIA configuration one adds D6 
branes on top of the NS5 branes. The strings between the D6 and the D4 branes play the role 
of flavored fundamental quarks of the SU(N) gauge theory [To]. T dualizing back to IIB the 
D6 branes get mapped into D7 branes that intersect the singular locus. The D7 branes in our 
configuration span the world volume coordinates of the D3 branes, the radial direction and wrap 
the three cycle of the deformed the conifold. The transverse coordinates are the coordinates of 
the S 2 part of the base. In terms of the transverse coordinates this means having one D7 brane 
at the north pole and another in the south pole of the S 2 . These "two" D7 branes are in fact 
one D7 brane since they smoothly connect at the origin of the radial direction. By writing down 
the 8d effective action of the D7 branes which includes the DBI term and the WZ term, we have 
shown that indeed these D7 brane configurations are solutions of the equations of motion. As a 
consistency check, the probe configuration satisfies the RR-flux cancellation condition. We will 
see that the probe configuration admits the isometry 577(2) x {7(1) which is the subgroup of the 
isometry group of the deformed conifold £{7(2) x SU(2). 

We then analyze the fluctuation modes around this classical D7 brane solution. The world- 
volume theory of the D7 probe contains U(Nf) gauge fields as well as (pseudo) scalars in the 
adjoint of the non-abelian flavor group. Kaluza-Klein(KK) reduction of the fluctuations around 
the compact three manifold of the 8d scalars yields 5d (pseudo) scalars whereas the vector fields 
yield vectors as well as pseudo scalars. We then determine the quadratic fluctuations which are 
solutions of the equations of motion of the 5d effective action. Both type of fluctuations have 
quantum numbers that are compatible with those of pseudo scalar and vector mesons of M = 1 
SQCD, namely they are in the adjoint of SU(Nf) and carry zero baryon number. To ensure that 
these modes are trivial under the £{7(2) x {7(1) global symmetry, we consider only a massless 5d 
vector fluctuation and the lowest scalar ones. 

A crucial point in the extraction of the mesonic spectra from the corresponding fluctuations 
is that these are normalizable with an appropriate regularity condition satisfied at the origin of 



2 



the deformed conifold. As a consequence, we show that the spectrum is discrete with a mass gap 
for the vector mesons. We also argue that the pseudo scalar meson spectrum exhibits the same 
behavior. It is shown that the system of mesons is stable (in the quadratic order approximation) 
since the kinetic term of the scalar fluctuations is positive and since there are no tachyonic masses. 

2 /3 

The mass scale of the mesons is that of the glueballs, namely, M meson ~ g € Na , . 

The paper is organized as follows. In section 2, the KS background is described in the 
alternative formulation of [TT]. The action of the D7 probe including the DBI part and the WZ 
is analyzed in section 3. We determine the solutions of the equations of motion that correspond 
to two D7 brane probes that merge smoothly into each other at the origin. Section 4 is devoted 
to the determination of the fluctuations around the probe configurations and the extraction 
of the mesonic spectrum. In subsection 4.1 the quadratic fluctuations are analyzed and then in 
subsections 4.2 and 4.3 the spectrum of the vector mesons and the scalar mesons is computed. We 
summarize and discuss some open questions in section 5. Some useful formulae are summarized 
in the appendix. 

2 The KS model and the geometry of a deformed conifold 

The KS background is based on adding fractional D3 branes into the deform conifold. The 
original KS metric jU] 



made use of the metric of a deformed conifold given in JU] • It turns out that for our purposes it 
is more convenient to use the formulation of ^T] since it admits a separation between the three 
cycle and two cycle of the deformed conifold. It is given by 



ds\ Q = h-^dxl + h^dxl 



(2.1) 




(sin uj 1 + cos oj 2 ) (d8) 



— (cos 9 cos (j) uj 1 — cos 9 sin <p uj 2 — sin 9 uj 3 ) (sin 9 d<p) 
H — K\t) sinh(r) dr 2 + (sin 6* cos 0a; 1 + sin # sin </> cu 2 + cos9u 3 ) 2 



(2.2) 



Here 



(sinh(2r) - 2r) 
2i/3 gink T 



1/3 



h(r) = (g s Ma') 2 2 2 / 3 e- 8 / 3 I(r), 



K(t) 



(2.3) 
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J(r) 



00 xcothx - 1 . 1/3 

ete 5 {smh{2x) — 2x) 

sinh x 



(2.4) 



For instance it is easy to verify in this formulation that for r = 0, the 6d metric reduced to 
J^(u; a ) 2 , giving the S 3 while the S 2 shrinks to zero. 

It is useful to rewrite the metric as follows 

dsj = R 2 [m 2 r x l 2 dxl + I^dil) , (2.5) 

where 



R 2 = 2 1/3 g s Ma', 



c 4/3 



(2.6) 



2 2 / 3 {g s Ma') 2 ' 

and e~ A / 3 ds\ = dx\. m is the mass scale of glueballs. For a detail of the glueball spectrum in the 
KS background, see [To] . 

The NS B-field is 

g s Ma> 



Bo 



(/ + k)g 3 A g A + V2kg 3 A d6 + V2kg A A sin0d(j) + 2k sinOdO A 



(2.7) 



Here 



f(r) 
k(r) 



r coth r — 1 

2 sinh t 
t coth r — 1 



(coshr — 1) , 
(coshr + 1) . 



2 sinh r 

For the definition of the one-forms g, see the appendix. The RR 2-form potential reads 

Ma' 
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- cos e'dip' A d<f)' + V2F(t) (g 3 A sin 6d<f> - g A A dO) 



Here 



Fir) 



sinh t — t 



(2.8) 



(2.9) 



(2.10) 



2 sinh r 

We will later need the explicit expression for C$ the RR six form. The corresponding RR 7-form 
field strength F 7 = dC 6 is defined by 



dF 7 = d*dC 2 - F 5 A dB 2 , 
where F5 is the self-dual 5-form field strength given by 

F5 — + *J r 5 , 

with 

T 5 = B 2 AdC 2 = ^-l(T)g 1 Ag 2 Ag 3 Ag i Ag 5 



(2.11) 



(2.12) 



a l{r) 



g s K 2 h 2 sinh r 
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(2.13) 



with a = A(g s Ma') 2 e- 8/3 and 

t cot n T — 1 

1(t) = f(l -F) + kF= — = (sinh 2r - 2r) . (2.14) 

4sinh t 

^From the equation of motion of C 2 it follows that dF 7 = 0. We thus obtain 

C 6 = —^—(f + k')h~ l d 4 x A dr (— sin 9 cos (f)uj l + sin# sin </>co> 2 — 2 cos9 d(fj 

- ^-(k- /)/i _1 d A x V2 (g 3 d6 + g 4 sin 6d(f)) + 2sin9d9d<j) . (2.15) 
where ' = 8/ dr. Finally we notice that the RR-4-form gauge potential is defined by 

F 5 = dC 4 + B 2 A dC 2 = B 2 A dC 2 + * {B 2 A dC 2 ) , (2.16) 

from which we obtain 

C A = 24/ ^ mR)4 u(r)d A x . (2.17) 



At) = T „J.\ 2 ■ (2-18) 



Here u(t) is defined by 

K 2 P sinhV 

Let us now review briefly some geometrical aspects of the deformed conifold. The deformed 
conifold is defined by 



e 2 



w±w 2 - W3W4 = — — . (2.19) 

For the explicit form of Wi, see the appendix. This can be regarded as a C* fibration over a 
two-dimensional complex plane spanned by w 3 ,w7 4 . The fibers get degenerate when W\W 2 = 0, 
giving a smooth singular locus in the base which is T-dual to an NS5-brane ^21- It turns out that 
the condition W\w 2 = can be solved as (1) 8 = 0, 6' = ir, (2) 9 = ir, 9' = n. Each corresponds 
to a cylinder that intersect with one another at a circle. To see this, note that for the two cases 
one finds 



wi =w 2 = 0, 

w-i = ■l—e^'"^ (± sinh - — cosh - 
V2 V 2 2 

w A = ^-e-W-V) f± sinh I + cosh I) . (2.20) 
\/2 V 2 2/ 

These denote two cylinders spanned by r, 0' — ip f that intersect with each other at a circle at 

r = 0. Recall that the circle is embedded in the S 3 and has a radius proportional to e. This 

circle corresponds to a diamond ^3]. In the limit e — >• where the deformed conifold reduces to 

a conifold, the radius goes to zero so that the smooth locus is split into two separate cones that 

intersect with each other at the tips of the cones. The T-dual of the conifold JH] gives us the 

IIA brane configuration that consists of two perpendicular NS5 branes, one at w\ = and the 

other at w 2 = 0. 
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Figure 1: Configuration of 7-branes on 2-sphere 



3 The D7 brane probes in the KS background 

Now let us consider D7 probes in the KS background that yield massive flavors, which were 
first discussed in [5]. For that purpose, it is useful to start from the T-dualized picture of the 
conifold limit e — > 0. The fractional D3-branes get mapped to D4-branes extending to the cycle. 
In this setup, one can add D6-branes on top of either of the two perpendicular NS5-branes to 
get Af = 1 SQCD with massless flavors [T7j, because the D4-branes divide the D6-branes into 
two pieces, each of which is responsible for chiral flavor groups. Upon T-dualizing the IIA brane 
configuration back, each type of D6-branes gets mapped to different D7-branes [18J: one is D7- 
branes which intersect with the singular locus W\ = and the other D7-branes intersect with 
w 2 = 0. 

Now turn on e ^ to have the deformed conifold. As seen before, the singular loci become 
a single smooth cylinder. Correspondingly, the two kinds of D 7-branes intersecting the singular 
loci become a smooth component of D7-branes that intersects with the smooth locus. More 
precisely, let us define D7(i)-branes staying at 6 = 0, the north pole of the S 2 and D7( 2 )-branes 
staying at 6 = ir, the south pole of the S 2 . See figure [TJ 

We show that these two seven branes intersect with each other at the S 3 at r = and form a 
smooth D7-brane: we first notice that the world volumes of the D7(i) and D7(2) are charactorized 



by 




(3 



= sin I e^'^') (sinh I T cosh I) , (3.1) 



V2 



and 



,(2) _ 
1,2 - 


e 


& 

cos — 
2 


( ( 2 ) _ 
3,4 - 




0' 


71 


sin — 
2 



iu$ = — sin - e^W -*)[- sinh - t cosh - ) . (3.2) 

It is easy to see that 

4\t = 0) = uf\r = 0) , d T 4\r = 0) = -9 r cf } (r = 0) . (3.3) 

These relations guarantee the smooth connection of the 7-branes at r = 0. 

Recall that in the deformed conifold the chiral flavor symmetry is broken to the diagonal 
subgroup U(Nf), as we are left with only a single type of D7-branes. This implies that turning 
on e amounts to turning on a vev for a Higgs field on the 7-branes that transforms as the 
bifundamental representation of the chiral flavor group. This vev plays the role of a mass term 
in the 4d gauge theory on the fractional D3 branes. 

Let us now analyze the D7-probe action in the KS background to verify that the configuration 
solves the equation of motion. The action consists of two parts 

S = Sdbi + Swz ) (3-4) 



where 



S 



= I d 8 a V- det (<f>*g + <j>*B 2 + 2na'F) , (3.5) 

S WZ = -fir I e^' F+B A J2 C P+ i 
J p 

= -fjL 7 J Q(2WF + 5) 3 C 2 + i(2WF + 5) 2 C4 + (2WF + 5)C 6 ^ . (3.6) 

Here <p*g,<f)*B2 are the pull-backs of g,B 2 . F is the U(l) field strength on the brane. fij = 
l/(27i) 7 a' 4 . The two transverse coordinates of a probe D7 are taken to be 9, (ft while upon 
taking the static gauge the world volume coordinates are given by X a = (X 1 ,9',(fi',ip f ), with 
X 1 = (x M , r). What we should show is that the D7(i) given by 9 = 0, any 0, and the D7( 2 ) given 
by 9 = it, any <fi solve the equation of motion. In order to check first the D7(i) brane, we assume 
9 = const, = 4>{X a ) and expand the action around = 0. It is verified that the coefficients 
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of the linear term in 9 depend linearly on some of the components of F. By requiring them to 
vanish, we find that the U(l) gauge potential on the probe should take the form 



A = A 3 (t)uj 3 . 

It is then not difficult to show that the action becomes 
jj 7 R s 



(3.7) 



S 



DBI 



9s 



J d 8 X^/-det(g {8) + R^B^) 



x 



16 ~r 

—At 



wz 



\i-jR 8 m A 
2*g, 



J d A xdTtuWuj 3 (l + 0{9 2 



1/2 



1 + 0(9 2 



2- 2 '\f + k)(f + k')r l + (2 5 /\f + k')!- 1 + 2 3 T ((f + jfe)u)) A 3 - 2 13 / 3 «' A 



(3. 



Here Aq 



and 



9(8)a/3 



7-1/2 

dX a dX p = gudX^X- 1 + —[K cosh r^u 1 ) 2 + {u 2 ) 2 ) + (Ksinh r)'(cu 3 ) 2 J , (3.9) 

(3.10) 



2d(0) 



/+* 1 A 2 



with 



Also 



gudX^X" 



jl/2 



m 4 



J = (if cosh r) 2 + x 2 



X 



(K sinh r)' c/r 2 . 
f + k 



(3.11) 



(3.12) 



2 4/3Jl/2 • 

Since there appear no linear terms in 9, we find that 8 = solves the equation of motion of 9. 
However A 3 = is not a solution because A 3 couples to a non-trivial current. We can easily see 
that the current is conserved. Recall that the contribution of the current from the DBI action 
is due to the non-trivial NS B-fields. 2 The equation of motion of A 3 is so complicated that no 
explicit form of the solution is found. The solution must be given to evaluate the tension of the 
probe 7-brane. We leave it as an open problem. However, the D7m probe configuration itself 
is allowed as the existence of the solution to the equation of motion should be guaranteed. It is 
interesting that the gauge configuration induces a D3-brane charge on the probe: 



F A F = 2 / rfrw'wV F t3 F 12 . 



(3.13) 



2 The backgrounds discussed in jS] have no non-trivial B-fields and hence the gauge field configurations on the 
probe branes are trivial. 
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s 2 

Figure 2: D7-probe configuration in a deformed conifold. The vectors denote the singular locus where the 
elliptic fibers get degenerate. The two shadowed surfaces that intersect with each vector are the 7-brane probes. 

Also by expanding the action around 6 = tt with any and the same form of the U(l) 
potential, we obtain the same form of the action as above. Hence the D7(2) probe also solves the 
equation of motion. 

Recall that (|3.9|) gives the metric of the 8d hypersurface that denotes the embedding of the 
D7(i) and D7(2) probes in KS. We will refer to the compact 3d submanifold as M3. M3 is topo- 
logically S 3 and admits the isometry group 577(2) x U(l), where SU(2) is the left transformation 
of T that leaves u a invariant and U(l) is the subgroup of the right SU(2) under which uj 1,2 
transform as a doublet while u 3 is invariant. One finds that the isometry group leaves _B 2 , C2, C4 
and A invariant as well. 

Here it is crucial to notice that the solution satisfies the consistency condition of the RR-flux 
cancellation. ^From the lOd point of view the two D7 branes carry opposite RR charges on S 2 
so that the total charge vanishes. See figure UJ 

To summarize, the probe D7-brane configuration looks like figure El 

4 Spectrum of mesons 

In this section, we study the fluctuation modes around the D7 probe configuration. Since these 
belong to the adjoint representation of U(Nf), so do the dual states in the QCD. We thus expect 
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to find the information about mesons made out of the fundamental quarks. Recall that no state in 
the adjoint of U(Nf) is charged under the baryon number C/ B (l), where U(Nf) = SU(Nf)xU B (l)- 
There exist two kinds of mesons we can learn about in the present context: pseudo-scalar and 
vector mesons. To see this, let us start with the D7 action defined on the probe configuration. 
We have two kinds of dynamical fields on it, 8d gauge potentials and 8d scalars that correspond 
to the fluctuation along the transverse directions. Upon KK reduction around M 3 , we obtain 
an effective 5d theory that consists of the infinite number of 5d vector and 5d scalar fields. We 
see that the 5d vectors are dual to vector mesons in the dual QCD, and the 5d scalars dual to 
pseudo-scalar mesons. Here we consider only the lowest-lying KK modes, namely, a massless 
5d vector and two 5d scalar fields, one coming from KK reduction of a 8d scalar by a scalar 
harmonics and the other from that of a 8d vector by a vector harmonics. The reason why we 
ignore the other massive modes is that these modes are irrelevant to the physics of QCD: since M3 
has the isometry SU(2) x U(l), the nontrivial harmonics carry non-zero spins and U(l) charges. 
However there is no counterpart of these quantum numbers in the dual QCD. 

In this section, we examine quadratic fluctuations of the D7 action with Nf = 1 around 
the probe configuration and derive the effective 5d theory that involves only the lowest-lying 
KK modes. As far as our concern is in quadratic fluctuations, we can concentrate on the case 
Nf = 1, since to this order any non-abelian interaction is irrelevant and consequently a field 
in the adjoint representation of U(Nf) reduces to Nj free fields. Based on this, we perform a 
numerical computation of the spectrum of vector and scalar mesons. In fact these are pseudo 
scalars rather than scalars as follows from a straightforward parity analysis. Recall that the 8d 
super Yang-Mills(SYM) theory on a D7 comes from dimensional reduction of lOd SYM. In ten 
dimensions, the parity transformation for the lOd vector field is defined as 



Upon dimensional reduction to 8d, we find that the 8d scalar and vector fields have parity odd. 
Since the scalar and vector harmonics on M 3 have parity even, the 5d scalars and vectors have 
parity odd. 

4.1 Quadratic perturbations around the probe configuration 

Here we consider the fluctuation around the D7(i)-branes in some detail. The analysis of the 
D7( 2 )-brane is then straightforward. In fact, one ends up with the same results in the two cases. 

As mentioned above, we assume 



VA(t,x 9 )V 



-1 



A(t, -Xg) . 



(4.1) 



9 = 9(X I ), = 0(X 7 ). 



(4.2) 
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The 8d scalars 9 and 4> reduce to the 5d scalars via a constant mode on M3. The constant mode 
is the normalizable zero mode of the scalar harmonics as M 3 is compact. The 8d vector A yields 
a 5d massless vector and a 5d scalar field. To see which vector harmonics is to be taken to obtain 
the 5d scalar, we first recall that the lowest-lying vector harmonics on S 3 transform under the 
isometry SU{2) x SU{2) as (3, 1) © (1?3). M 3 at hand is topologically S 3 and admits as the 
isometry the subgroup SU(2) x £7(1). The lowest-lying vector harmonics then splits into some 
of the vector harmonics on M3 that belong to an irreducible representation of SU(2) x £7(1), one 
of them being singlet. In fact, it is found that the singlet vector harmonics is given by uj 3 . We 
thus find that the the 5d scalar corresponds to the fluctuation of A% in ()3.7|) around the solution 
of the system ()3.8|) . However, since the solution is not known, one can not obtain the explicit 
form of the action that governs the fluctuation. As we will see below, up to quadratic order this 
fluctuation mode decouples from the other 5d scalar and vector modes. In this paper, we focus 
on only the two modes to compute the meson spectrum. 

It turns out that up to quadratic terms in the fluctuations, the DBI action takes the form 

/i 7 i? 8 m 4 



S 



DBI 



J ^Xlu'lu 2 ^ (iTsinhryj 1 / 2 



x 



i + L g iJ dl e dj e + l 



2W 



g IJ g KL F IK F JL 



+ 



2W 



i/V 3 ^3^ 3+ ™^ 



Al - 



2W 8 



x 



R 2 I l l 2 J 



A, 



Here 



K cosh t 
4I l / 2 J 



IK 2 smh 2 r + ^(rcothr - 1) 



(4.3) 



(4.4) 



a,b = 1,2,3 are the indices of the orthonormal frame of S 3 , see the appendix for detail. u a b is 
defined by the 3d part of the tensor + R~ 2 B^ and given by 

/ i^coshr —x \ 



jl/2 



Uab 



V 



—X 

+x K cosh r 
(fTsinhr)' J 



(4.5) 



On the other hand, we find that the fluctuations from the WZ term up to quadratic order are 
given by S'wz in ()3.8|) . This shows that 9 and Aj have no mixing terms. As mentioned before, 
A3 can not be diagonalized because one has to expand this around a non-trivial configuration 
that solves the action (|3.8|) . 

It is interesting to notice that the kinetic term of 9 is always positive. This implies the 
stability of the probe configuration at hand. 
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4.2 Vector mesons 



Let us first study the 5d vector part to compute the vector meson spectrum. Part of the analysis 
given here is parallel to JH] that aimed at investigating a localization of a bulk gauge field at a 
brane world. 

The action of the 5d gauge potential is 

S v ~ m 4 J d 5 X (K sinh r)' \fl g IJ g KL F IK F JL . (4.6) 

Here we neglect the overall numerical factor because this is irrelevant for the purpose of computing 
the mass spectrum of mesons. The equation of motion reads 

dj ({K sinh t)' g IJ g KL F IK } = 0. (4.7) 

For the component L — r, this becomes 

= f3,(3A-M,). (4.8) 

It is useful to work in the gauge A T = 0. Then the above equation gives us the Gauss law 
constraint. We next decompose the rest components of the gauge potential in terms of the 
complete set of functions Xnij)'- 

A M (x,r) = X)4 n) ^)Xn(r). (4.9) 

n 

X satisfy a second-order differential equation we will present below. Note that we are interested 
in only normalizable solutions. Substituting the decomposition into the Gauss law constraint, 
we obtain 

= J2^ u d,A^(x)d TXn (r). (4.10) 

n 

As we will see later, the constant mode is not normalizable so that we find 

f^=0. (4.11) 
This implies that Ajf 1 are Proca fields. Substituting (|4.9J) into the action, we obtain 
S v = Jd 4 xdr (K sinh r)' \fl 

E [lrrirF$FW X n{T)x m (T) + {K l^ hr)/ rf v AfA^ d TXn (r) d rXm (r)_ ■ 

n,m ^ 1 

(4.12) 



12 



Now we define \n as the solution of the differential equation 



yp) \{K sinh r)'/(r 
with the normalization condition given by 



' d r [ -n? • ^777 ^ rXn ) = A nX n(r), (4.13) 



(4.14) 



o 



Here 



^7 = (Ksinhr)'I{T) y/J^rj . (4.15) 
It then follows that the action becomes 

S V = J2 [ d ^ x V^V^F^Fj^ + ^K^A^A^ . (4.16) 

Thus the Proca fields A satisfy the on-shell condition with the mass square given by 

M 2 n = 4A n m 2 . (4.17) 

We regard this as the mass square of vector mesons of QCD. 

The differential equation ()4.13|) allows two independent solutions: one is normalizable and 

the other non-normalizable. We are interested in the normalizable solution here. It turns out 
that the normalizable solution should behave as 

X(r) = e-*7(T), (4.18) 

with the boundary behavior of / 

f(r -> oo) = const. (4.19) 
As a check, we notice that the einbein ^py behaves as 

^7~(r -> oo) ~ r 3/2 . (4.20) 



This guarantees the normalization condition (|4.14|) . 

It is easy to verify that / obeys the differential equation 

f" (r) + A(r)f'(r) + (B(r) + XC(r))f(r) = 0, (4.21) 

with 

1 4 
A = - d T log J - - , 

B = - l -d T \ogJ+^ , 

C = (KfanhT)'I(r) . (4.22) 
13 



We will solve this differential equation numerically following the procedure discussed in [20] to 
find the glueball spectrum of QCD. We first find out the asymptotic behavior of the solution 
at r ^> 1 for a generic A. Using this data as an input, the solution can be found numerically. 
By imposing a regularity condition at r = to be discussed in a moment, only solutions with 
appropriate values of A are allowed. In particular, we will see that the spectrum of vector mesons 
is discrete with a mass gap. 

In order to obtain the asymptotic solution, we notice that for large r 

A{t) = *n(r)e- 2n \ B{t) = Ur)e- 2nT , C{r) = e"T £ c n {r)e~ 2n \ (4.23) 

n=0,l,- ra=0,l,- n=0,l,- 



where 



2 50 8r 

100 , 16r 



0, h 



9 ' 9 
4r 2 32r 3 



c = r, ci = — , c 2 = — , • • • (4.24) 



Expanding /(r) as follows 



/(r)= E ( 4 - 25 ) 

71=0,1,- 

it follows from (|4.21|) that the coefficients f n obey the recursion relation 

/n - y/n + -g-/n + ^ ° m ( fn-3m ~ 3 ( n ~ 3m ) /n-3m ] + &m/n-3t 

m=0 L ^ ' 



+ A E Cmfn-Sm-l = • (4.26) 



m=0 



By setting / = 1, the solution is given by 

f - 1 f - 9AT f - 27AV f - 25 8U3r3 f4 27) 

/o-l, /i g ' /2 "-6^' /3 "12"^02T' ( 7) 

Now let us discuss what is the regulatory condition to be imposed at r = 0. As seen before, 
the probe D7 brane consists of the two pieces, D7(i) and D7( 2 ). So far we have examined the 5d 
massless vector potential only on one of the two 7-branes. The true solution is given by interpo- 
lating smoothly between the two solutions each of which solves (|4.13jl . We denote X™ ( r )j Xn (t) 
by the solutions on D7m and D7(2), respectively. Since both obey the same differential equation 
with the same asymptotic behavior, the two solutions are related as 

x£ ) (t) = ±xZHt). (4.28) 
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This shows that the regularity condition is given by 

drXn(r = 0) = or Xn (r = 0) = 0. (4.29) 

It follows from the regularity condition of vanishing of the derivative of \n that the allowed values 
of the eigenvalue A are 

A = 0.50, 1.64, 3.66, ■ •• (4.30) 
which give us the mass spectrum of vector mesons 

M 2 

— - = 2.00, 6.56, 14.6, (4.31) 
m 1 

and from the regularity condition of vanishing \ n that the values of A are 

A = 0.86, 2.50, 4.98, ••• (4.32) 

which corresponds to vector mesons masses 

M 2 

— - = 3.44, 10.0, 19.9, (4.33) 

m 2 

4.3 Pseudo scalar mesons 

In this subsection, we compute the scalar meson spectrum. As discussed in the previous subsec- 
tion, there exist two scalar modes in the 5d effective action that correspond to the spectrum of 
scalar mesons that carries no SU(2) x U(l) quantum numbers: one is 9 in (J4.3)) . and the other 
is A 3 . 

Let us first consider 9. As shown in ()4.3j) . the fluctuation on both of the probe 7-branes is 
governed by 

5 SC ~ m 2 J d 4 xdr(K sinh t)'Vj ^g^drf dj9 . (4.34) 
We first decompose 9{X T ) in terms of the complete set of appropriate functions £ n (r). 

6(x»,T) = J20 {n) (x)UT), (4-35) 

n 

By defining £ as the solution of 

-^'{W^rfl 8 ^)'^' <436) 
with the normalization condition given by 

I dTy/p^ n (r)^ m (r) = 5 nm , (4.37) 
J o 
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where 

y/p=(KanhT)'y/IJL, (4.38) 

the action becomes 

S sc = 1 I d 4 x [rTd»0 {n) dvO (n) + 4a ri mV n) # (n) ] . (4.39) 

n 

Thus we obtain the scalar mesons with the mass square 

Ml = Aa n m 2 . (4.40) 
In order to have the normalizable solutions, we see that that £ should behave as 

£(r) = g(r) , (4.41) 

with 

g(r — > cx)) = const. (4.42) 
As a check, we note that the einbein ^/p behaves as 

v^(r -> cx)) ~r 3/2 e 2r/3 . (4.43) 

We find that g obeys the differential equation 

g"(r) + L>(r)</(r) + (£(r) + aC(r))g(r) = 0, (4.44) 

with 

D = d T log ' ^ 



(Ksinhr)'// 3 ' 

£ = (4.45) 

As before, we first solve the asymptotic behavior of g for a generic a. For that, we need the 
asymptotic behavior of the coefficients D, E: 

D = E e ~ 2nrrf «> E e ~ 2nre «> ( 4 - 46 ) 

n=0,l,--- n=0,l,--- 



where 



w 4 w - 8r w - 32r ' 

«o — — g ; _ ~3~ ' ~ 3~ ' ' " 

32r 128r2 ,„ „N 

e = ei = — — , e 2 = — , • • • (4.47) 
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_ 2nT 

(Tie 3 . 



Expanding g(r) as 

n=0,l,- 

it follows from ()4.44|) that the coefficients g n obey the recursion relation 



(4.48) 



. An . 4n 2 ^ 
9 n ~ Y 9n + ~g~ 9n + 

m=0 



d m ( g'n-3m ~ 3 ( n ~ 3m ) 9n-3m 



n—3m 



3m— 1 



(4.49) 



m=0 



By setting ^ — 1> the solution is given by 



#0 = 1, 9i 



3ar 



92 



27a 2 r 2 



9?> 



It 81ctV 



3^-3 



(4.50) 



4 ' J ~ 128 ' 15 2560 

Using this, we solve ()4.44j) numerically. As before, we have to impose the regularity condition at 
r = 0: 

f ( T = 0) = , or <9 T £(r = 0) = . (4.51) 

We obtain from this 

a = 1.77, 3.91, 6.90, ■■• (4.52) 
which give us the mass spectrum of scalar mesons 



M 2 

— - = 7.08, 15.6, 27.6, 



(4.53) 



Unlike the vector mesons, it turns out that for the pseudo scalar ones the same values of a 
are found for both types of the regularity condition. This may be related to the fact that around 
r = due to the shrinking of the S 2 the fluctuations and hence £ vanish anyhow. 

The computation of the meson spectrum associated with A% is left as an open question. 



5 Discussion 

In this paper, we have discussed adding D7-brane probes to the KS background for the purpose 
of getting a SUGRA description of M = 1 SQCD with flavors. The point here is that for 
Nf N c the backreaction of the D7-branes to the KS can be suppressed. To find the probe 
configuration, the geometrical data of the deformed conifold and its T-dual brane picture in 
IIA were useful. Based on this observation, we discussed the configuration and verified that it 
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solves the equation of motion of the probe action. As a consistency, we argued that the probe 
configuration satisfies the RR-flux cancellation condition. Using this result, we next computed 
the scalar and vector meson spectrum by examining the normalizable fluctuations around the 
probe with an appropriate regularity condition at r = 0. We discussed that the mass spectrum is 
charactorized by the single mass scale m, being equal to the glueball mass, and has a mass gap. 
An important open problem here is to compute the scalar meson spectrum that comes from A 3 . 

There still remain some issues to be explored. One is to check whether the probe configuration 
preserves M = 1 supersymmetry by examining K-symmetry on the probe brane, although that is 
plausible from the T-dualized IIA brane picture. It would be interesting to generalize the probe 
configuration such that one is allowed to have one parameter family of the solutions. For instance, 
D7 and D5 probe configurations in AdS§ x S 5 were found [S] that depend on one parameter, being 
identified with the mass of flavors. It is also nice to compute Wilson lines via SUGRA to see 
how the screening of a pair of heavy quarks in terms of dynamical quarks occurs, as was done for 
AdS§ x S 5 with D7-probes |Sj. Baryons could be realized as D3 brane probe which is wrapping 
the M 3 and is connected with N strings to the D7 brane probe. 

In this paper, we have been working in the probe approximation, which is justified for Nf <C 
N c . In order to have a SUGRA configuration for any Nf, we need to find a fully localized D7- 
brane configuration in the KS. To achieve this sounds rather difficult, however. In fact, only a 
few example of fully localized brane solution are known [21]. Instead of working in the KS, it is 
useful to consider the Penrose limit of the KS jllj . where the original KS background gets much 
simplified and the string spectrum on it is simple enough to work out in light-cone gauge. One 
may expect to obtain a fully-localized D7-brane solution in the resultant background. It would 
be interesting to study the string theory on that background. 

As we have seen, the D7 probes yield massive flavors whose mass is of order of the glueball 
mass. Our model is not similar to a realistic QCD, unfortunately. In fact, we are not allowed 
to take naively the limit e — > to obtain massless flavors, because this limit gives rise to a 
singularity in the KS and therefore the SUGRA approximation is not valid any more. It is quite 
interesting to find a probe configuration in an well-defined SUGRA background that provides us 
with a more realistic model with light quarks with chiral flavor symmetry. 
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A Formulae 

Here we summarize some useful relations. 
Consider 

T = e^ V3 e^ Vl e^ V3 , S = e^ a3 e^ 6ai . (A.l) 

T defines an S 3 as a Hopf fibration over S 2 spanned by 9 f , <fi' and S defines an S 2 . u a are defined 
as the left-invariant one- forms of the S 3 : 

T ] dT=-u a a a . (A.2) 



2 



Explicitly 



w 1 = cosib'd9' + sin 9' sin if)' d(j)' , 
w 2 = sinip'dO' — sin 9' cos ip' d(f) , 

w 3 = dip' + cos 9' sinip'dcf)'. (A.3) 



The left-invariant one-forms obey 



du a = -e abc u b u c . (A.4) 



2 

The volume form is 

w 1 A w 2 A w 3 = sin & deb 1 A d& A dib' . (A.5) 

The dual vectors are 

_ ,, d sin-0' d n > ■ 1/ 9 

o cos ib 1 d , . d 

* = smi ' m - w + 0!ii ' W ' 

* = w ■ (A ' 6) 

In order to find the relation between the basis used in [5] and that in the present paper, one 
first carries out a change of basis for the 1-forms uj % [11J, 

yfifo x -yfit°* + tfo* = S^ a a a S, (A.7) 
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which gives 



g 5 = sin 9 cos (f) U 1 — sin 9 sin <fi u 2 + cos 9 uj 3 , 
1 



-9 



V2 
1 



(sin^u 1 + cos (f)uj 2 ), 



g 3 = —= (cos 9 cos 4> cu 1 — cos 9 sin (fiuj 2 — sin 9 to 3 ). 
V2 



After a little algebraic work one finds 



9 5 = 9 5 

e-^(g l +ig 2 ) = (g 3 - v^sintf #) + % (g 4 + yftdO) 



(A.8) 



(A.9) 



We next denote Wi that defines the deformed conifold W1W2 — W3W4 = — y in terms of S and 
T. Define 

L = TS, R = a 3 Sa 3 , (A.10) 

from which, we define 



W = LW e lfi=( Wl W3 ), W e = -^ei< 



This gives the deformed conifold: 



One can find that 



detW = . 

2 



(A.ll) 



(A.12) 



Wl 2 



U>3A 



e 

it 
V2 



sinh - cos 9 cos -e^'+V-') T l sinh 1 sin sin ^ e ±! W-V) ± cosh 1 cos °l e ±W+V) 
2 2 2 2 2 2 

sinh - cos 9 sin ^W'V) _ sinh - sin cos ^ e ±*(W+VO ^ cosh 1 sin ^ e ±§(<^') 
2 2 2 2 2 2 



(A.13) 
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